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Abstract —In this paper, a novel method for sensor node local¬ 
ization under mixed line-of-slght/non-line-of-sight (LOS/NLOS) 
conditions based on second order cone programming (SOCP) 
is presented. SOCP methods have, hitherto, not been utilized 
in the node localization under mixed LOS/NLOS conditions. 
Unlike semldefinlte programming (SDP) formulation, SOCP is 
computationally efficient for resource constrained ad-hoc sensor 
network. The proposed method can work seamlessly in mixed 
LOS/NLOS conditions. The robustness of the method is due to 
the fair utilization of all measurements obtained under LOS 
and NLOS conditions. The computational complexity of this 
method is quadratic in the number of nearest neighbours of the 
unknown node. Extensive simulations and real field deployments 
are used to evaluate the performance of the proposed method. 
The experimental results of the proposed method is reasonably 
better when compared to similar methods in literature. 

Index Terms —Semidefinite programming. Second order cone 
programming. Localization, Ad-hoc sensor networks, LOS and 
NLOS conditions 

I. Introduction 

Localization is a crucial aspect in vehicular sensor networks, 
where link instability among sensor nodes in NLOS conditions 
poses a significant challenge. The problem of localization 
using optimization techniques is generally non-linear and non- 
convex in nature. Hence, a closed form solution to the local¬ 
ization problem is not always achievable. The problem can 
be relaxed using semidefinite programming (SDP) or second 
order cone programming (SOCP) into a convex formulation. 
The sensor node localization problem using SDP and SOCP 
has been extensively dealt with in [l]-[5] and [1], [6], [7] 
respectively. A relation between SOCP and SDP formulation 
is analysed in [6]. It is shown here that SDP provides a 
tighter relaxation for the sensor node localization problem 
when compared to SOCP. In the presence of noise, both SOCP 
and SDP relaxation provide suboptimal solutions [6]. There is 
a trade-off between SDP and SOCP in terms of computational 
complexity and localization error [5]. Notably, the efficacy of 
solving the localization problem using SOCP relaxation lies 
in least computational complexity than the counterpart SDP 
[6]. The localization problem using SOCP can be solved in 
polynomial time. Hence, for tiny sensor nodes with resource 
constraints, it is imperative to choose a computationally power 
efficient method. 

A classical SOCP formulation for distributed sensor node 
localization is proposed in [5], [7]. The extension of distributed 
SOCP (D-SOCP) method [7] to account for NLOS conditions 
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increases the localization inaccuracy of sensor node. SOCP 
method for sensor network localization with anchor uncer¬ 
tainty is suggested in [6]. However, though this method con¬ 
siders anchor uncertainty, it fails to capture NLOS conditions. 
This often occur in practice. A non-cooperative and coopera¬ 
tive localization based on SOCP and SDP relaxation without 
any explicit assumption of NLOS is discussed in [1]. SDP 
approach for sensor network localization with noisy distance 
measurements is presented in [4]. However, [4] considers only 
measurement noise in distance measurement. SDP algorithm 
for sensor node localization [3] is limited to uncertainties in 
anchor positions and propagation speed. This method does 
not consider the case of NLOS positive bias. The problem of 
localization under NLOS condition is addressed in [2]. This 
method is based on SDP relaxation, which is computationally 
expensive compared to SOCP. 

In this work, we propose a SOCP based relaxation method 
for sensor node localization in Mixed LOS/NLOS conditions 
(MLN-SOCP). Non-linear least square [3] and approximate 
maximum likelihood [2] objective functions are exploited 
in the context of SDP formulation. Subsequently, a NLOS 
localization problem is formulated in terms of SOCP [1], 
[6], [7] for lower computational complexity. The proposed 
method is shown to work in polynomial time. Other major 
advantages of proposed method are its ability to work in mixed 
LOS/NLOS conditions [8]-[10], and robustness to noise under 
random node deployment scenario. The proposed method 
provides fairness to all distance measurements obtained from 
both LOS and NLOS. In a large sensor network, there is 
a scarcity of anchors and obtaining the measurements may 
become expensive, if NLOS measurements are discarded. 

The notations used in this paper are as follows. The bold 
faced letter (upper or lower case) represents the matrix. (.)^ 
denotes the transpose of a matrix. The Euclidean distance 
between two one-dimensional vectors is represented by I 2 
norm, ||.||. The cardinality of a set is denoted by |.|c. The 
radio communication range of a node is R in Nd x Nd 
network. The rest of the paper is organized as follows. Section 
II describes the SOCP formulation for sensor node localization 
problem. Performance evaluation is demonstrated in Section 
III, and a brief conclusion is presented in Section IV. Solution 
to the localization problem in SEDUMI form is listed in the 
Appendix. 

11. SOCP Method for Node Localization in Mixed 
LOS/NLOS Conditions 

In this Section, problem formulation for the sensor node 
localization is described first. A SOCP relaxation is described 
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next. Subsequently, computational complexity and conver¬ 
gence analysis of the proposed method are presented. 

A. Problem Formulation 

The sensor network is modelled as undirected topology, Q = 
(H,£,n7). The set of nodes including anchors is represented 
by H = {1, 2,..., rriu, ■ ■ ■ ,Po}- The first nodes are un¬ 
known sensor nodes, H^, while next pa — mu are anchors 
or known sensor nodes. Ho. The coordinates of the unknown 
node and anchor are denoted by € R'^ (d = 2 or 3) and at 
€ R'^, V r,t respectively. The set of edges, £ = £i U £„ among 
nodes is assigned a weight, '[Ur,t S w called link quality . £i 
and £n represent the set of LOS and NLOS edges respectively. 
The channel is assumed to be invariant in both the forward and 
reverse direction i.e. Wr,t = ‘^t,r ■ 

The line-of-sight (LOS) distance between receiving node 
and emitting node (anchor or known node) is given by 

d\ t = ||xr - ajll -f Wr,t, (r,t) e£i,r < t, (1) 

where Wr^t ^ is the measurement error, which 

follows a zero-mean Gaussian distribution with variance cr^j. 

The non-line-of-sight (NLOS) distance between receiving 
node and anchor is corrupted by a large positive error, t 
due to obstruction. 


and rju respectively. In order to obtain the location of unknown 
sensor nodes, Equationj^needs to be minimized over the entire 
network. This optimization problem is non-linear in nature and 
it is difficult to find a closed loop solution. Finding the global 
optima is also not easy, because the problem is non-convex in 
nature. 


B. Solution using SOCP Relaxation 

To solve this optimization problem, second order cone pro¬ 
gramming (SOCP) relaxation is used, which is computation¬ 
ally efficient for sensor node localization. Defining auxiliary 
variables {qr,t} and {Zr,t}, the problem in Equation |5] can be 
reformulated using relaxed conic constraints as 


minimize 


+ ^r,t) 

{r,t)GS;r<t 


subject to 



- atll - 4.t| < {r,t) € £i 

I 11 d^ t I ^ (f, f) C £n 


(6) 

Introducing an epigraph variable V [ 6 ], [13] and auxiliary vari¬ 
ables U and {yr,t], the equivalent convex epigraph problem 
is written as 


= ||xr - at|| -b Or,t + Wr,t, {r,t) e£n,r < t, ( 2 ) 


where Or-,t follows a exponential distribution with mean pLr,t 
and variance ^ [2], [11], [12]. A similar model in the context 
of semidefinite programming is described in [2]. The NLOS 
distance is simplified by subtracting the mean, assuming 
known NLOS state [2]. Consider 7 ^^ = p'^ t + a^ t, 

t^r.i = - Pr,t = ||Xr - atll +w'^t, {r,t) €£n,r < t, (3) 

where w',. ^ is approximately modelled as Gaussian random 
variable, Af{0,^rt) [ 2 ]- For large error, this approximation 
may suffer from the long tail problem. However, it is valid 
for the most of the scenarios as described in [2]. The joint 
distribution of LOS and NLOS measurement distances is given 
by 


= n [dr.t'PuiKdlt -\K - atW) 

{r,t)G£;r<t 

|LOS} + (1 - gr,t)Pw'{{d'ft - ||x, - a^lDlNLOS}], 

( 4 ) 

where Vui and Vw' are the probability density function of 
the error Wr^t and w'^t- The probability of LOS link is 
denoted by Following [2], the objective function using 
the approximate maximum likelihood method can be recast as 


minimize 

{r,t)GS;r<t 


2 9r,t |||^r ^t\\ ^r,t 


l^r.t 


1 I 12 

H ^ (1 9r, t) 11 11 t 

Ir.t 


( 5 ) 


The variance of LOS link is = rif{d\.tY, while variance 
of NLOS link is 7 ^^ = {rff -b Standard deviation 

per unit length of noise and NLOS error are represented by rji 


minimize V 


subject to 


< V 


\/^ jyr,t - dl tj < qr,t, (r,t) € £l (7) 

V 


' 1 - gr,t 


\yr,t i I ^ ^ 


|Xr 11 ^ yr,t 


where U = Fourth constraint is induced from 

the relaxation of equality constraint into an inequality one. 
Localization of the sensor node is now reduced to solving 
Equation [T] which represents a SOCP relaxation of original 
non-convex optimization problem. The solution of Equation |7] 
is obtained by numerical optimization technique like SEDUMI 
solver [14]. The solutiorQis enumerated in Appendix. 


C. Computational Complexity Analysis 

Let 2 *^ unknown sensor node can be localized using pi 
number of neighbouring anchors within its vicinity, which are 
contained in set Pi. The computational complexity [1], [15] 
of the proposed method is given by 



2 = 1 2 = 1 


where e and M denote the total number of equality constraints 
and iteration complexity respectively, kj and N represent the 
dimension of second order cone and total number of 

’Solution to the problem in SEDUMI form is listed in the Appendix 










3 


oi 0.16 



Fig. 1; Figure illustrating the (a) Cramer-Rao lower bound and 
(b) Root mean squared error plot for the proposed method in 
a 40 m X 40 m network; for p = 0.0025, gr^t = 0.7, |H|c = 
1600, R = Nd, rji = 0.1, rjn = 0.06. 

second order cone constraints respectively. The total number of 
conic constraints is 3pi +1, V f as can be noted from Equation 
|2l The number of equality constraint for the SOCP formulation 
in Equation |2] is zero. The dimension of first to fourth set of 
constraints are {2pi + 1), 2, 2 and 3 respectively. Therefore, 
asymptotic computational complexity for the proposed MLN- 
SOCP method is reduced to 0{M\Pi\'^), which is quadratic in 
cardinality of Pi. Computational complexity can be reduced 
by constraining \Pi\c to localize an unknown sensor node. 

III. Performance Evaluation 

In this Section, experimental conditions for real field de¬ 
ployment is presented first. Subsequently, experimental results 
for sensor node localization in mixed LOS/NLOS are dis¬ 
cussed. 

A. Experimental Conditions 

An experimental ad-hoc sensor network is deployed using 
both National Instrument (NI) WSN node and Crossbow 
motes to validate the effectiveness of the proposed algorithm. 
For a larger dimensions of the network, NI WSN - 3202, 
3212 and gateway NI 9792 are used. For a smaller network 
dimensions. Crossbow MTS310 sensor board, MIB520 USB 
mote interface board and XM2110 IRIS board are used. The 
nodes communicate among themselves through IEEE 802.15.4 
protocol. Sensor nodes are randomly deployed inside the 
network assuming a uniform distribution. 

B. Experimental Results 

In this Section, Cramer-Rao lower bound analysis is illus¬ 
trated first. Subsequently, localization error analysis followed 
by probabilistic error analysis of location estimation are pre¬ 
sented. 

1) Cramer-Rao Lower Bound Analysis: In order to assess 
minimum variance of the localization error, Cramer-Rao lower 
bound [16], [17] analysis plot is illustrated in Eigure [T] Eor 
sake of illustration, anchors are placed at the boundary of 
the network, while nodes are deployed inside the convex hull 
of networks. The minimum CRLB on localization error is 
found to be 0.073, which is at the center of the surface. 
The corresponding root-mean-square error (RMSE) plot using 


estimated distance from the proposed algorithm is shown in 
Figure [Itb). The minimum value attained for the proposed 
method, MLN-SOCP is 1.4719, while for D-SOCP method 
[7] is 3.2682. 

2 ) Localization Error Analysis: In this Section, variation of 
localization error with various network parameters namely net¬ 
work dimension, Nd, radio communication range, R, fraction 
of anchors, p and probability of LOS links, are illustrated. 
Localization error is defined as the Euclidean distance between 
actual and estimated node location. 

(a) Comparison of Localization Error with Varying Network 
Dimension and Radio Communication Range: If the anchors 
are linearly separable over all nodes in the networks, localiza¬ 
tion error [18] is given by 

E||x,-x,|| =0(iv|i?§|HU^) (9) 

where network is considered to be of size Nd x Nd- The 
communication range of the node is R, while |H|a represents 
the cardinality of set of anchors. E represents the expectation 
operator. As the network dimension increases, localization 
error increases as shown in Table U The proposed method, 
MLN-SOCP performs significantly better than the D-SOCP 
for all values of radio communication range, R. The increase 
in the localization error with network dimension can be at- 
most linear [18]. However, if the ratio of network dimension 


TABLE I; Comparison of localization error (m) for various 
Nd, R, Tji and pn, for p = 0.3, gr,t = 0.7, |H|c = 100. 


Methods 

Radio 

Range 

rji = 0.10. rin = 0.06 

Nd = 40m 

Nd = 40m 

Nd = 80m 

m = 0.2 
Tin = 0.15 

VI = 0-3 

Vn = 0.25 

D-SOCP 

V2N^ 

4.30 

10.05 

6.99 

8.40 

Nd 

3.99 

8.34 

6.97 

7.90 

MLN-SOCP 

V2Nd 

3.84 

8.24 

4.40 

6.52 

Nd 

2.83 

4.27 

4.39 

6.33 


to the number of anchors is kept constant, then it is likely 
that localization error is unaffected. On the other hand, if the 
number of anchors are held constant, there is an increase of 
localization error. This is due to the fact that the node may not 
able to hear sufficient number of anchors for the localization. 
Increase in radio communication range enables a node to hear 
many anchors. But, incorporating many anchors increases the 
possibility of inclusion of erroneous distances between anchor 
and node. This results in an increase of the localization error 
at the rate of C9(i?3) [18]. There is a marginal increase in 
localization error with standard deviation of noise and NLOS. 

(b) Comparison of Localization Error with Varying LOS 
Probability and Node Density: As the node density increases, 
the localization error reduces as expected in Table HI! However, 
localization error decreases till a certain |H|c. Increasing the 
node density implies an increase in number of anchors by a 
proportion ’p’. Inclusion of more anchors leads to the increase 
of erroneous distance between unknown node and anchor. 
This results in reduced localization accuracy after, |H|c = 200. 
There is a marginal decrease in localization error on increasing 
the LOS probability links, as expected. 
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TABLE II: Comparison of localization error (m) for varying 
|H|c and for p = 0.3, R = Nd = 40m, pi = 0.1, rjn = 
0.06. 



1 '^1 c 

Methods 

grit 

50 

100 

150 

200 

250 

300 

D-SOCP 

0.95 

4.05 

3.70 

3.66 

3.54 

3.91 

4.01 

0.7 

4.36 

3.99 

3.78 

3.65 

4.09 

4.49 

0.4 

5.68 

5.52 

5.16 

5.11 

5.27 

5.73 

0.1 

5.79 

5.56 

5.45 

5.25 

5.31 

5.88 

MLN-SOCP 

0.95 

2.64 

2.55 

2.48 

2.39 

2.44 

2.53 

0.7 

2.85 

2.83 

2.54 

2.50 

2.74 

2.77 

0.4 

2.87 

2.85 

2.60 

2.59 

2.81 

2.82 

0.1 

2.91 

2.90 

2.69 

2.68 

2.84 

2.86 


3) Probabilistic Error Analysis of Location Estimation: 
The cumulative distribution function (CDF) represents the 
statistical distribution of localization error. Figure |2] shows the 
empirical CDF computed for various values of p. Since the 
proposed MLN-SOCP method uses the maximum likelihood 
model, the MLN-SOCP performs significantly better than the 
D-SOCP method for higher values of probability of anchors. 
This aspect can be verified from Figure |2l since the CDF curve 
rapidly attains a probability of unity. However, for lower values 
of p = 0.1, MLN-SOCP method performs reasonably better 
than D-SOCP. 



Fig. 2: Figure illustrating the CDF plot for various p; for 
= 0.7, |H|c = 100, R = Nd = 40m, pi = 0.1, pn = 0.06. 


IV. Conclusion 

In this paper, second order cone programming for sensor 
node localization in mixed LOS/NLOS is proposed. The 
novelty of the proposed method lies in the development of 
NLOS model using SOCP relaxation. The experimental results 
on sensor node localization in mixed LOS/NLOS illustrate that 
the node can be localized with a high degree of reliability 
even for large number of NLOS links using the maximum 
likelihood model. The results also indicate that the proposed 
method can be effectively used for node tracking in vehicular 
ad-hoc networks. New methods that can utilize SOCP to detect 
and localize malicious anchors are currently being explored. 

V. Appendix 

Let unknown sensor node can be localized using pi 
number of neighbouring anchors within its vicinity. 

S,r=[xr Vr Xl X2 ... Xp, 2/pJ [2pi-r2) X1 


The objective function of Equation |7] can be written in SE- 
DUMI form as 


where 


X = [H^ C^r,t Zr.t 

yr,t V f 

(11) 




’^r,t = \_^r,l 2 r,2 

• ^’'■P*]lxpi 

(12) 

yr,t = \yr,i yr,2 

y^’Pi]lXpi 



The objective function in the SEDUMI is of the form 

B'^x = V, where Bt = -B = - [ 0 (i_ 5 p.+ 2 ) l]^ (13) 

The total number of constraints are 3pi + 1. The dimension of 
first to fourth set of constraints are 2pi -f 1, 2, 2 and 3 respec¬ 
tively. Comparing with standard form of conic constraint, 

||A"^x- f c|| < b^x-f d (14) 

For first constraint in Equation |2l 

A^x = U, Cl = 0, b^x = V and di = 0 (15) 


Therefore, 


Ai = [Ei E 2 


Efc 


E„ 


where 


Efc = [Efc^ Efe_ 2 ] 

Efe,l = [ 0 (i_ 2 pi-|- 2 ) Ifc.t 0 (l, 2 pi-|-l) 
Efc, 2 = [ 0 (l, 3 pi-l- 2 ) Ifc.t 0 (l,Pi-l-l)] 
br = [0(i,5p.+2) 1 ] 


(16) 

(17) 

(18) 


with lfc,t is row vector of length pi with a 1 corresponding 
to the index k and O’s elsewhere. Similarly, for second set of 
constraints 

iT 


A 2 ,fc — [0(i_4p^_|_2) 


Lfc,t 


0 ]- 


b 2 ,fc = 


Vi 


--dk,t [0(i, 


2pi-|-2) ife.i 


s/ 9 k,t 

C2,fc = [0 , C2 = [C2,l 

A2=[[b2,l A2 ,i] ... [b2.. 


0 


(i, 2 pi-ri) 

^2,Pi] 


(19) 


Pi A2,pi]] 


On similar line, matrices for third set of constraints can be 
written. For fourth set of constraints, 

1 T 


A4,fc = 


1 0 
0 1 


(l, 5 pi-l-l) 

(l, 5 pi-l-l). 


b 4 ,fc = [0(4pi-r2,i) lfc,t 0] , C4,fc = [0 -af ] 

Similarly, A 4 and C 4 can be written like A 2 and C 2 . 

A = -[[bi Ai] A 2 ... A 4 ] 

C = [Oi,p, C2 C3 04 ]^ 


( 20 ) 

( 21 ) 

( 22 ) 


Cartesian product of all cones which is the dimension of 
constraints 


Q — [2pi + 1 2.1(1 2pi) 8.1(i,p.)] (23) 

Hence, Equation |7] can be written in terms of A, Bt , C and 

Q 
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